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Аннотация 


Введение. Компьютерное моделирование позволяет инженерам принимать обоснованные проектные решения за 
счет точной оценки тепловых характеристик объектов проектирования. Актуальным направлением научных ис- 
следований и разработок является реализация технологии цифровых двойников в процессе проектирования тех- 
нических объектов. Для этого необходимо разрабатывать компьютерные модели, точность которых соответ- 
ствует требованиям, предъявляемым к цифровым двойникам. Однако в научной литературе недостаточно ши- 
роко представлены результаты исследований, направленных на реализацию технологии цифровых двойников в 
процессе проектирования. В основном рассматриваются общие вопросы, связанные с применением цифровых 
двойников в различных отраслях промышленности. Поэтому целью данного исследования явилась разработка 
цифровой модели и сравнительный анализ точности расчетов тепловых характеристик объекта проектирования. 
Материалы и методы. В качестве основного инструмента для проведения исследования выступает предложен- 
ная авторами методика разработки компьютерной модели тепловых характеристик для реализации технологии 
цифровых двойников. Численное решение реализовано путем построения тепловой модели для расчета темпера- 
турного поля на основе метода конечных элементов в системе инженерного анализа «Ansys» от компании «Ansys 
Inc» (США). Для аналитического решения применяется разработанная на основе метода пространства состояний 
компьютерная модель тепловых характеристик, реализованная в модуле «Ansys Twin Builder». Модель простран- 
ства состояний приводится в соответствие с поведением исходной тепловой модели путем приближения переда- 
точной функции к пошаговому отклику тепловой нагрузки с применением метода векторной аппроксимации во 
временной области. Верификация построенной аналитической модели выполнялась в системе инженерных рас- 
четов «MATLAB» от компании «The MathWorks» (США). Исследования проводились для станка модели 400V 
производства предприятия ООО «НПО «Станкостроение» г. Стерлитамак (Россия). 

Результаты исследования. Разработана цифровая модель, позволяющая с высокой точностью выполнить расчет теп- 
ловых характеристик объекта проектирования. Результаты сравнительного анализа показывают высокую степень со- 
ответствия значений тепловых характеристик, полученных с помощью предложенной цифровой модели, результатам 
численного моделирования. Максимальная погрешность расчета тепловых характеристик не превышает 0,1 °С. 
Обсуждение и заключение. Компьютерное моделирование, сочетающее численные методы расчета и научный 
подход, основанный на технологии цифровых двойников, позволяют получить результат максимально прибли- 
женный к результатам экспериментов. Предложенная в исследовании цифровая модель является эффективным 
решением, поскольку позволяет выполнить расчеты для оценки тепловых характеристик в режиме реального 
времени, что является одним из важнейших требований при реализации технологии цифровых двойников. 


Ключевые слова: цифровой двойник, сложный технический объект, компьютерное моделирование, 
автоматизированное проектирование, температурное поле, тепловые характеристики 
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Abstract 

Introduction. Computer modeling allows engineers to make valid design decisions by accurately assessing the thermal 
characteristics of design objects. The implementation of digital twin technology in the process of designing technical 
facilities 1$ the current direction of scientific research and development. To do this, it is necessary to develop computer 
models whose accuracy meets the requirements for digital twins. However, the scientific literature does not widely present 
the results of research aimed at implementing digital twin technology in the design process. The general issues related to 
the use of digital twins in various industries are mainly considered. Therefore, the objective of this study was the 
development of a digital model and a comparative analysis of the accuracy of calculations of thermal characteristics of 
the design object. 

Materials and Methods. The main tool for conducting the research was the methodology proposed by the authors for 
developing a computer model of thermal characteristics for the implementation of digital twin technology. The numerical 
solution was implemented through constructing a thermal model for calculating the temperature field based on the finite 
element method in the ANSYS engineering analysis system from ANSYS, Inc. (USA). For the analytical solution, a 
computer model of thermal characteristics developed on the basis of the state-space method, implemented in the ANSYS 
Twin Builder module, was used. The state-space model was matched to the behavior ofthe original thermal model through 
approximating the transfer function to the stepwise response of the thermal load using the time domain vector 
approximation method. Verification of the constructed analytical model was carried out in the engineering calculation 
system MATLAB from the MathWorks company (USA). The research was carried out for a 400V machine model 
manufactured by NPO "Stankostroenie" LLC, Sterlitamak (Russia). 

Results. The developed digital model makes it possible to calculate the thermal characteristics of the design object with 
high accuracy. The results of the comparative analysis showed a high degree of correspondence between the values of 
thermal characteristics obtained using the proposed digital model and the results of numerical simulation. The maximum 
error in calculating thermal characteristics did not exceed 0.1°С. 

Discussion and Conclusion. Computer modeling that combines numerical calculation methods and a scientific approach based 
on digital twin technology, provides obtaining the result as close as possible to the results of experiments. The digital model 
proposed in the study is an effective solution, since it provides performing calculations to evaluate thermal characteristics in real 
time, which is one of the most important requirements for the implementation of digital twin technology. 


Keywords: digital twin, complex technical object, computer modeling, computer-aided design, temperature field, 
thermal characteristics 
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Введение. Компьютерное моделирование традиционно является эффективным инструментом для решения 
тепловых проблем на ранней стадии проектирования сложных технических объектов. Однако решение про- 
блемы требует точной оценки тепловых характеристик объекта проектирования для снижения негативных эф- 
фектов, вызванных повышением температуры [1]. При этом одним из эффективных инструментов для предва- 
рительной оценки тепловых характеристик является имитационное моделирование в системе инженерного 
анализа на основе передовых цифровых решений и разработок. Например, в работе [2] авторов Jianying Xiao и 
Kaiguo Fan представлен разработанный цифровой двойник для определения тепловых характеристик. Принцип 
работы двойника заключается в моделировании теплового поведения объекта путем отображения и коррекции 
тепловых граничных условий. Результаты эксперимента показали высокую точность модели (более 95 96), что 
существенно для повышения точности моделирования тепловых характеристик и тепловой оптимизации. По- 
этому актуальным направлением научных исследований и разработок в области моделирования является при- 
менение систем искусственного интеллекта [3] и цифровых двойников [4]. Например, в работе [3] авторов 
Наогап У1 и др. для повышения точности анализа тепловых характеристик предложена интерактивная модель 
коррекции тепловых граничных условий на основе нейронной сети. Результаты экспериментов показывают, 
что точность расчета температурного поля превышает 98 %, а точность прогнозирования тепловой деформа- 
ции — 96 %, что эффективно повышает точность моделирования. Кроме того, в работе [5] авторов Kypra- 
нова Н.В. и др. отмечено, что цифровые двойники часто применяются с целью совершенствования физических 
прототипов сложных технических объектов, поскольку позволяют не только обеспечить информационное со- 
провождение процесса проектирования, но и принимать эффективные проектные решения на основе наработок 
в области искусственного интеллекта. 

Одной из характерных особенностей технологии цифровых двойников является то, что для имитационного мо- 
делирования часто применяются компьютерные модели пониженного порядка [6]. Поэтому разработка компьютер- 
ных моделей является одним из базовых условий в реализации технологии цифровых двойников [7]. Понижение 
порядка моделей — эффективный и понятный с математической точки зрения подход для преодоления ограничений 
времени выполнения многомерных имитационных моделей. Например, в работе [8] авторов Мирзаев Д.А. и др. 
предложена простая аналитическая модель тепловых полей для разработки цифровых двойников процесса про- 
мышленной дуговой сварки. Авторы Бордачев Е.В. и Лапшин В.П. в своей работе [9] представили результаты ма- 
тематического моделирования температуры в зоне контакта инструмента и изделия при токарной обработке метал- 
лов. Такой подход позволяет получить точную оценку тепловых характеристик, соответствующую результатами 
численных экспериментов в режиме реального времени. В работе [10] авторы Schróder С. и Matthias У. представили 
модель пониженного порядка и предложили новую процедуру балансировки модели, основанную на преобразова- 
нии сдвига состояния. В заключении представили результаты сравнительного анализа с результатами из источников 
литературы, полученных с помощью серии численных экспериментов. 

Применение линейной и инвариантной по времени компьютерной модели пониженного порядка позво- 
ляет обеспечить быстрое моделирование при сохранении высокой точности расчетов [11]. При разработке 
компьютерной модели выполняется аппроксимация [12] передаточной функции для приближения модели 
пространства состояний к пошаговому отклику исходной тепловой модели [13]. Поскольку пошаговый от- 
клик тепловой нагрузки получен из исходной тепловой модели, то цифровая модель должна предоставить 
те же значения тепловых характеристик. 

Однако несмотря на то, что в последнее время наблюдается рост интереса к цифровому двойнику, в научной 
литературе недостаточно широко представлены результаты исследований, связанных с реализацией цифровых 
решений в процессе проектирования технических объектов. На основе систематического обзора литературы и 
тематического анализа публикаций по цифровым двойникам выявлен один из ключевых пробелов в знаниях, 
связанный с разработкой математического, программного и методического обеспечения высокоточных компью- 
терных моделей в рамках реализации технологии цифровых двойников. 

В этой связи целью исследования являлась разработка компьютерной модели пониженного порядка и анализ 
эффективности её применения в составе цифровой модели для точного расчета тепловых характеристик сложных 
технических объектов проектирования. 

Для достижения поставленной цели требовалось построить тепловую модель и выполнить расчет темпера- 
турного поля объекта проектирования, сгенерировать независимые пошаговые отклики тепловой нагрузки с по- 
мощью разработанного программного сценария [14], реализовать цифровую модель для расчета тепловых харак- 
теристик, определить погрешность расчетов, полученных с помощью численного и аналитического решений и 
провести сравнительный анализ полученных результатов моделирования. 
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Материалы и методы. Построение температурного поля объекта моделирования выполнялось для однород- 
ного изотропного тела на основе уравнения нестационарной теплопроводности: 


OT[0t = a? AT +9, с, (1) 


где T — температура (°С); t — время (c) a -4/A/(p-c) — коэффициент температуропроводности (м/с); 
X— коэффициент теплопроводности (Вт/м-°С); с — удельная теплоемкость (Дж/кг:°С); p— плотность материала (кг/м?); 
А — оператор Лапласа; д, — объемная мощность тепловыделения (Br/M?). 

Тепловой поток в процессе теплопередачи принимался равным количеству теплоты, переносимой через про- 
извольную поверхность площадью 5 в единицу времени f, что выражается следующим уравнением: 


бы | f аза, Q) 


где О — тепловой поток (Br); qn — плотность теплового потока (Br/w?); 5 — площадь поверхности (M°). 

Плотность теплового потока при теплоотдаче определялась по формуле: 

qe -9-S-(Ts -Т.), (3) 

где a — коэффициент теплоотдачи (Вт/(м?-°С)); Ts — температура поверхности (°С); To — температура окру- 
жающей среды (°С). 

В этой связи, для расчета температурного поля объекта моделирования по формуле (1), задавались тепловые 
потоки (2) и (3), которые определяли количество теплоты, проходящее через поверхность в единицу времени. 

В качестве объекта моделирования был выбран компонент (рис. 1) металлорежущего станка модели «400V» 
производства предприятия ООО «НПО «Станкостроение» (Россия, г. Стерлитамак) в виде сверлильной головки. 


Рис. 1. Геометрическая модель объекта в «Ansys Design Modeler»: 
1 — корпус; 2 — электродвигатель; 3 — шпиндельный узел; О — тепловые потоки; 
qv — объемная мощность тепловыделений; qn — плотности тепловых потоков; 
a — коэффициент теплоотдачи; Ti-T4 — температурные датчики 


Поскольку для расчета температурного поля учитывалось количество теплоты, выделяемое в основном элек- 
тродвигателем (электромагнитные потери) и шпиндельным узлом (механические потери на трение), то в качестве 
внутренних источников тепла принимались электродвигатель qyi и шпиндельный узел 4,2, вблизи которых зада- 
вались соответствующие тепловые потоки О! и Q». Плотности тепловых потоков qni, 92 назначались для поверх- 
ностей, расположенных вблизи электродвигателя и шпиндельного узла соответственно. 

Конвекция задавалась коэффициентом теплоотдачи а с учетом условий теплообмена (свободная конвекция B 
воздухе). Поскольку станок контактирует с газообразной средой (воздух), то количество теплоты, отданное 
нагретой поверхностью окружающей среде в единицу времени f, прямо пропорционально разности температуры 
поверхности Ts и среды To в зависимости от 5 площади теплоотдающей поверхности (3). 

При построении тепловой модели объекта (твердого тела), состоящего из однородного материала (конструк- 
ционная сталь) с постоянными теплофизическими свойствами и наличием внутренних источников тепла, назна- 
чались следующие начальные и граничные условия: 

— начальные условия учитывали фиксацию постоянной температуры по всей поверхности объекта моделиро- 
вания (t= 0: T= То = const); 

— граничные условия второго рода задавались тепловыми потоками электродвигателя (Q1), шпиндельного узла (Q2), 
плотностью теплового потока (9.1) от электродвигателя к передней стенке корпуса и (4,2) к внутренним поверхностям; 

— граничные условия третьего рода задавались коэффициентом теплоотдачи (Q) для поверхностей, располо- 
женных внутри корпуса сверлильной головки; 

— граничные условия четвертого рода для контактных соединений поверхностей учитывали идеальный теп- 
ловой контакт и отсутствие термического сопротивления: 


Ta zum +4,. (4) 


+0 ? 
On |0 On |, 
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Тепловые потоки (Qi, Q2), плотности тепловых потоков (qni, 42), как и объемная мощность тепловыделения 
(qvi, 9), назначались в соответствии с известными рекомендациями по металлорежущим станкам [1]. Граничные 
условия, как и тепловые потоки, задавались для наружных поверхностей. Поэтому учитывалось, что взаимосвязи 
тепловых полей присутствовали только между внешними поверхностями. 

Дифференциальное уравнение (1), совместно с начальными и граничными условиями второго (2), третьего (3) 
и четвертого (4) рода, представляет собой математическую формулировку поставленной задачи. Поставленная 
задача решалась с помощью методов численного и аналитического моделирования. 

Численное решение выполнялось на основе метода конечных элементов в системе инженерного анализа 
«Ansys», которая разрабатывается американской компанией «Ansys Inc» (США) и поставляется фирмой «МЦД», 
авторизованным дистрибьютором «Ansys» в России. Геометрическая модель (рис. 1) объекта импортировалась в 
проект «Ansys Workbench» с добавлением блока анализа переходного температурного поля «Transient Thermal». 
При численном решении выполнялась калибровка параметров моделирования и граничных условий тепловой 
модели (таблица 1) для приближения модельных значений температуры к экспериментальным данным. 


Таблица 1 
Граничные условия (параметры) тепловой модели 


" Плотности тепловых Коэффициент 
Мощность тепловыделений| Тепловые потоки ов а 
Параметр SUOTA 


qvi, Вт/м? 42, Вт/м? Qi, Вт Q2, Вт qni, Вт/м? qw, Вт/м? а, Br/(w?-?C) 


Значение 6 500 1 000 28 15 32 18 15 


В модуле «Ansys Mechanical» назначались начальные (начальная температура 790-24 °С) и граничные (таб- 
лица 1) условия для построенной сеточной модели, выполнялось построение температурного поля (рис. 2). При 


этом коэффициент теплопроводности принимался равным 4=60,5 Вт/(м °С) и назначался таковым для конструк- 
ционной стали. 


Тепловая модель объекта включала два контактных соединения для электродвигателя и гильзы шпинделя с 


корпусом сверлильной головки и содержала 7 тепловых граничных условий. Построенная сеточная модель со- 
стояла из 16 309 элементов и 58 527 узлов. 


Общее время моделирования 21 600 секунд (6 часов) разбивалось на интервалы (1 час) и шаги Af = 360 секунд 


(6 минут), всего № = 60 шагов, в пределах которых параметры тепловой модели (граничные условия) принима- 
лись как постоянные и независящие от времени. 


Моде! (В4) 
Hy Geometry 
#- „8 Materials 
Ey Coordinate Systems 
9-и) Connections 
#-- „2 Mesh 
=-- „М Transient Thermal (B5) 
vT- Initial Temperature 
Mit Analysis Settings 
5-29 Simplorer 
/®. Internal Heat бепегабоп1 
/@. Internal Heat Generation2 
у. Heat Flow1 
sa Heat Flow2 
85 HeatFiuxi 
85 Heat Flux2 
Vf Convection 
S- Solution (B6) 
m- Solution Information 
s Temperature 
Bg «09 Simplorer 
8 ТетрРгоБе1 
8 TempProbe2 
T TempProbe3 
„® TempProbe4 


Рис. 2. Температурное поле объекта B «Ansys Mechanical» 


Расчет температурного поля требовался для генерации независимого пошагового отклика тепловой нагрузки, 
содержащего информацию об изменении температуры по времени (тепловые характеристики) для каждого пара- 
метра (таблица 1) тепловой модели в отдельности: 

УЕ (t) i=1,k, j=1,m, (5) 


где k — количество температурных датчиков; т — количество параметров тепловой модели. 
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Генерация пошагового отклика выполнялась в модуле «Ansys Mechanical» с применением расширения 
«Application Customization Toolkit (ACT)», которое поддерживает внедрение пользовательских сценариев. Это 
позволило разработать программный сценарий на языке программирования «Python» с целью автоматической 
генерации специального набора файлов, содержащих независимые пошаговые отклики [14]. 

Аналитическое решение выполнялось на основе метода пространства состояний путем построения модели 
тепловых характеристик по формуле: 

ї = Ах+ Ви 


, 6 
y = Cx + Du (6) 


где х = 0х / ôt — производная вектора х состояний по времени f; y и и — векторы выходных и входных данных 


соответственно; А, В, C, D — матрицы постоянных коэффициентов. 


В уравнении (6), вектор x-—(xi3»...xw) содержит переменные состояния, вектор входа 
u = (qui, 9, Оп, Q2, Чт, qm, Toi, Too, ..., Ток)! значения параметров тепловой модели (граничные и начальные усло- 
вия), вектор выхода y = (Ti, 75, ..., Т) значения тепловых характеристик. 


Модель передаточной функции выражается следующим уравнением: 


т Е 
„(= У, Hu()uy. io Lk, (7) 
где Н — матричная комплексная передаточная функция. 

Матричная передаточная функция Н получена путем применения к формуле (6) преобразования Лапласа, что 
выражается следующим уравнением: 

H(s)-C(sI- A4) Вв+р, (8) 
где 5 — комплексная переменная Лапласа; / — единичная диагональная матрица. 

Передаточная функция (8) отражает зависимость преобразования Лапласа выходной переменной 
Y(s) = H(s)U(s) от преобразования Лапласа входной переменной U(s) модели (6) при нулевых начальных 
условиях х(0) = xo = 0. При этом размерность матрицы передаточной функции H зависит от величины выхода 
к= 4 и размерности входной величины т = 10, что соответствует размерности исходного пошагового OT- 
клика (5). Поэтому модель (6) приводилась в соответствие с поведением исходной тепловой модели путем 
приближения ее передаточной функции (7) к пошаговому отклику (5) с применением метода векторной ап- 
проксимации [15]. 

Для построения матриц коэффициентов модели (6) выполнялся обратный переход от передаточной функ- 
ции (8) к модели в пространстве состояний. В этом случае, передаточная функция (8) принимает вид уравнения, 
в знаменателе которого содержится характеристический полином степени / = 4 (порядок системы), а в числителе 
полином степени Z = [- 1: 

z І 
G(s) = Е Р(5) = bg ОЗЕ Q(s) = dg 2232077 (9) 
где a u b коэффициенты полиномов Q(s) и P(S) соответственно. 

Корни полиномов Q(s) и P(s) представляют собой полюса и нули передаточной функции (8) соответственно. 
Метод неопределённых множителей применялся к уравнению (9) для разложения каждого элемента матрицы Н 
на элементарные дроби. Обозначив полюса характеристического полинома через p;, получаем уравнение следу- 


G()-D« Y". А 


i 
, 
5 — р; 


ющего вида: 


(10) 


2. Р(5) 
где R; E e 


Paur матрицы Ri обозначался через F; и выполнялось ее разложение на произведение двух матриц с полным 


— матрица размерностью (k х т); 4 — количество полюсов. 


рангом столбца и строки соответственно: 
kxr; pnxm =: 
R; 2 C?" Bj". rank(R;)— rj. (11) 


Матрицы модели (6) являются диагональными, размерностью А", В", (Cen. [)" (п = 56) и содержат эле- 
менты, которые получены непосредственно из коэффициентов передаточной функции (10). 
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Матрицы A системы и В управления содержат следующие коэффициенты: 


ра 0 ^а 0] bj 0 .. 0 
0 i 0 0 b e s. O 
P Pr „В , (12) 
0 0 рз 0 0 b, bj, 
10 0 0 Pan | [0 O0 ... bym 
Матрицы С выхода и D прямой связи содержат следующие коэффициенты: 
Ca e Сы 0 sis 0 о... 0 0 .. 1.. 0 0 
C- alt m ме е нн 3 
0 es 0 (кл) C(k—1)g 0 ... 0 0 ue Oa uu 1 0 ( ) 
о... 0 0 ses 0 Ck; e Cpg 0. .. 0 .. 0 1 


Такой подход позволил выполнить переход от модели B форме передаточной функции (8) к модели B пространстве 
состояний (6). Алгоритмы перехода также реализованы в системе инженерных расчетов «MATLAB», коммерческой 
компании «The MathWorks» (США), в виде специальных функции «ss2tf()» и «tf2ss()». 

Для получения значений тепловых характеристик решалась задача Коши для системы обыкновенных диффе- 
ренциальных уравнений, поскольку в формуле (6) переменная X является производной от вектора состояний тем- 
пературного поля по времени t. Решение системы уравнений (6) было получено с помощью метода Рунге-Кутты 
четвертого порядка. 

Верификация построенной аналитической модели выполнялась путем проведения серии вычислительных экс- 
периментов с применением, разработанной в системе «MATLAB» прикладной программы, которая включала pe- 
ализацию метода Рунге-Кутты четвертого порядка (рис. 3). Вычислительные эксперименты проводились на пер- 
сональном компьютере (процессор AMD Ryzen 5 56000 with Radeon Graphics 2.30 GHz, оперативная память 
16,0 ГБ, тип системы 64-разрядная операционная система Windows 10 Pro версия 21H2), характеристики которого 
являются базовыми для современной вычислительной техники. 

Построенная аналитическая модель и применение метода Рунге-Кутты четвертого порядка для решения си- 
стемы дифференциальных уравнений позволило с высокой точностью вычислить значения тепловых характери- 
стик (рис. 3). Поскольку максимальные значения погрешности, то есть разность между модельными значениями 
тепловых характеристик, полученных с помощью численного (FE-Model) и аналитического (LTI-Model) реше- 
ний, не превысили 0,72 °С на всем интервале моделирования. 

Линейная и инвариантная по времени (Linear and Time-Invariant, LTI) компьютерная модель пониженного порядка 
(Reduced Order Model, ROM) разрабатывалась в модуле «Ansys Twin Builder» системы инженерного анализа «Ansys» 
на основе температурного поля объекта и сгенерированного в модуле «Ansys Mechanical» пошагового отклика. ccc 
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Рис. 3. Графики тепловых характеристик в системе «MATLAB» 
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Реализация цифровой модели тепловых характеристик на основе температурного поля с применением век- 
торной аппроксимации заключается в последовательном выполнении семи основных этапов: 

Этап 1. Импорт геометрической модели объекта и построение тепловой модели в системе инженерного анализа; 

Этап 2. Расчет температурного поля объекта на основе построенной тепловой модели; 

Этап 3. Генерация независимого пошагового отклика по результатам численного моделирования тепловых 
характеристик объекта; 

Этап 4. Применение алгоритма векторной аппроксимации для получения полюсов и нулей передаточной 
функции модели пространства состояний; 

Этап 5. Построение модели пространства состояний по известной передаточной функции; 

Этап 6. Разработка компьютерной модели тепловых характеристик; 

Этап 7. Реализация цифровой модели объекта. 

Таким образом предложенная цифровая модель, содержащая компьютерную модель для точной оценки теп- 
ловых характеристик сложных технических объектов проектирования, была получена путем последовательного 
выполнения всех представленных выше этапов. 

Результаты исследования. Разработанная компьютерная модель (Thermal 56400У SML) содержит 6 вхо- 
дов и 4 выхода (рис. 4), что позволяет выявить взаимосвязь между параметрами тепловой модели и значениями 
тепловых характеристик. В качестве входных данных компьютерной модели приняты объемная мощность теп- 
ловыделений (qui, 4), тепловые потоки (Qi, Q2) и плотности тепловых потоков (41, 4,2). В качестве выходных 
данных выступают тепловые характеристики (T;—74). Компьютерная модель входит в состав цифровой MO- 
дели (рис. 4 а) тепловых характеристик, реализованной в модуле «Ansys Twin Builder». 

Значения параметров цифровой модели, представленные в табличном виде (рис. 4 в), подаются на вход ком- 
пьютерной модели (рис. 4 a) с помощью компонентов «STEP». Графический модуль демонстрирует значения 
тепловых характеристик (рис. 4 6), полученных на выходе компьютерной модели. 


Thermal SG400V SML 
RC 
34,0- 
T, 1 
34,01 
T; ] 
T, 2904 
Т, ] 
о аррар 
0 5 000 10 000 15 000 Lc 
a) 6) 
1 2 3 4 5 6 
Time, c 3 600 7 200 10 800 14 400 18 000 21 600 
qn2 18,0 18,0 18,0 18,0 18,0 18,0 
qvl 6 500,0 6 500,0 6 500,0 6 500,0 6 500,0 6 500,0 
01 28,0 28,0 28,0 28,0 28,0 28,0 
02 15,0 15,0 15,0 15,0 15,0 15,0 
Qv2 1 000,0 1 000,0 1 000,0 1 000,0 1 000,0 1 000,0 
qnl 32,0 32,0 32,0 32,0 32,0 32,0 


Рис. 4. Реализация цифровой модели в системе «Апзуз Twin Builder»: 
а — компьютерная модель тепловых характеристик; Ó — модуль графического представления результатов; 
в — модуль контроля входных данных компьютерной модели 


При разработке компьютерной модели задавались пределы размерности от 2 до 4 порядка, значения целевой 
ошибки £ = 5x10? и допуск для нулевого порядка £o = 2x10?. Остальные параметры устанавливались автомати- 
чески, поскольку метод векторной аппроксимации максимально автоматизирован в сравнении с остальными ме- 
тодами, которые поддерживаются в модуле «Ansys Twin Builder». 

В процессе разработки компьютерной модели модуль «Ansys Twin Builder» автоматически генерирует специ- 


альную матрицу ошибок аппроксимации М; = || у;()-у;(0 ||; во временной области, каждый элемент которой 
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отражает разность между значениями тепловых характеристик пошагового отклика (5) и передаточной функ- 
ции (7) модели в пространстве состояний, что выражается следующим уравнением: 
0,43 4,04 1,53 4,97 2,22 1,82 
- 021 2,22 2,54 1,07 1,16 1,95 103 (14) 
2,66 1,63 0,39 0,26 3,43 1,82 
1,46 3,21 0,06 0,84 2,95 3,50 

В данном случае максимальная относительная ошибка не превышает значение £ = 4,97х103. При этом все 
остальные значения ошибок оказались меньше указанного 5 = 5x10? предела. Нулевое значение ошибки в мат- 
рице (14) означает, что вход был проигнорирован в связи с очень малым вкладом. 

Оценка точности расчета тепловых характеристик с применением предложенной цифровой модели выполня- 
лась путем сравнительного анализа результатов, полученных с помощью численного и аналитического решений. 
Сравнительный анализ проводился по критерию максимальной погрешности. Максимальная погрешность AT, 
то есть разность между значениями тепловых характеристик, полученных с помощью численного и аналитиче- 
ского решений для всех температурных датчиков, в каждый момент времени рассчитывается по формуле: 

AT nax = тахАТ,, (15) 


тах 
j-ln 


где АТ; = |Т,/- T;4 — погрешность; Т; ги 7; и — значения температуры конечно-элементной и цифровой моделей 


соответственно ( = 1, т); т — количество температурных датчиков. 

Для оценки точности цифровой модели выполнялся расчет погрешностей, результаты которого представлены 
в виде поверхности (рис. 5 а) и линейного графика (рис. 5 6) максимальной погрешности по времени. Поверх- 
ность (рис. 5 a) представляет собой расчетные значения погрешности AT; для каждого температурного датчика 
(ат, аТ2, dT3, dT4) в отдельности и в каждый момент времени на всем интервале моделирования. 
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Рис. 5. Погрешности расчета тепловых характеристик: а — погрешность для каждого температурного датчика; 
Ő — максимальная погрешность для всех температурных датчиков 
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График (рис. 5 6) показывает, что для выбранного интервала времени максимальная погрешность AT, не 
превышает 0,1 °С. Погрешность расчета тепловых характеристик для каждого температурного датчика в отдель- 
ности (рис. 5 а) также не превышает заданных пределов. 

Обсуждение и заключение. Полученные результаты вычислительных экспериментов и проведенный срав- 
нительный анализ подтверждают эффективность применения предложенной цифровой модели, которая позво- 
ляет с высокой точностью AT. = 0,052 °С выполнить расчет тепловых характеристик, для дальнейшего прове- 
дения процедур анализа и идентификации тепловой модели. 

На температурное поле объекта проектирования влияет множество факторов, что осложняет определение номи- 
нальных значений тепловых граничных условий. Для решения данной проблемы в исследовании сначала анализи- 
руются ключевые технологии, задействованные в реализации цифрового двойника сложного технического объекта, 
затем строятся тепловая модель и компьютерная модель пониженного порядка LTI ROM переходного температур- 
ного поля. Разработанная компьютерная модель применяется в составе цифровой модели, которая позволяет полу- 
чить точную оценку тепловых характеристик объекта проектирования и тем самым повысить эффективность вы- 
полнения проектных процедур в процессе автоматизированного проектирования сложных технических объектов. 

Точность и эффективность вычислений компьютерной модели пониженного порядка и исходной тепловой 
модели полного порядка оценивается путем сравнительного анализа результатов моделирования. Результаты вы- 
числительных экспериментов показывают, что с точки зрения точности расчетов компьютерные модели пони- 
женного порядка и конечно-элементные модели полного порядка в основном сопоставимы по точности, а макси- 
мальная погрешность расчетов находится в пределах допустимого диапазона и не превышает 0,1 °С. 

Полученные в ходе вычислительных экспериментов результаты не противоречат результатам, представленным в 
источниках научной литературы по схожей тематике и позволяют сделать вывод о том, что применение предложенной 
цифровой модели эффективно для оценки тепловых характеристик сложных технических объектов в режиме реаль- 
ного времени, что является одним из важнейших условий для реализации технологии цифровых двойников. 

Однако изменение температурного поля сложного технического объекта по-прежнему зависит от многих фак- 
торов. Поэтому в дальнейших исследованиях необходимо разработать компьютерные модели тепловых дефор- 
маций и ввести эффективные алгоритмы оптимизации на основе искусственного интеллекта, чтобы обеспечить 
надежность результатов моделирования, полученных с помощью цифровых двойников. 
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